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,  20.  Abstract  (continued) 

are  described  and  compared  to  data  from  the  recovered  experimental  model.  The 
tunnel  closure  predictions  and  the  predicted  locations  of  block  motion  are 
in  reasonable  agreement  with  some  experimental  data. 

Other  aspects  of  the  study  included: 

•  development  and  application  of  an  analytical  technique  for 
predicting  joint  activation  conditions  (i.e.,  block  motions) 
for  specified  joint  orientations  and  locations  near  the 
tunnel  surface, 

•  calculations  and  comparisons  of  tunnel  response  to  dynamic 
and  quasi-static  stress  loads,  and 

•  investigation  of  the  effects  of  joint  friction  on  tunnel 
damage. 
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SECTION  1 


INTRODUCTION  AND  SUMMARY 


1.1  BACKGROUND 

All  rock  masses  contain  large  numbers  of  planar  discontinu¬ 
ities  or  joints.  These  joints  often  react  to  an  applied  stress 
quite  differently  than  the  surrounding  material  and  there  is  both 
experimental  (Ref.  1)  and  theoretical  (Refs.  2,  3)  evidence  to 
suggest  that  sliding  motion  across  these  joints  should  be  one  of 
the  important  considerations  involved  in  designing  underground 
structures  which  are  to  survive  nuclear  and/or  conventional 
explosions . 

In  the  recent  DIABLO  HAWK  underground  nuclear  test,  SRI 
International  fielded  several  laboratory-scale  structures  in 
jointed  rock  simulant  (Ref.  4).  Figure  1  shows  the  general  model 
designs.  The  jointed  rock  models  are  18M  and  30"  in  diameter, 
with  3 "  anu  5"  diameter  tunnels,  respectively.  The  SRI  small- 
scale  jointed  structures  experiment  provides  an  opportunity  for 
examining  the  influence  of  highly  idealized  jointed  arrays  on 
tunnel  closure  and  damage  characteristics. 

There  are  both  continuum  and  explicit  approaches  for  pre¬ 
dicting  tunnel  closure  and  damage  in  jointed  media.  In  the 
continuum  approach,  the  overall  mechanical  properties  of  the 
field  are  synthesized  from  the  properties  of  the  intact  rock  and 
the  properties  and  arrangement  of  the  joints;  for  example,  stiff¬ 
nesses  and  strengths  used  in  the  continuum  model  are  lower  than 
the  intact  rock  properties.  In  the  explicit  approach,  the  prop¬ 
erties  and  arrangement  of  individual  blocks  and  joints  are 
individually  and  explicitly  modeled  and  used  to  calculate  tunnel 
response  to  loading.  Both  of  these  treatments  have  important 
.limitations:  The  rationale  for  choosing  a  set  of  degraded 


COARSE  JOINT  SPACING  IN  OUTER  ANNULUS  UNIFORM  JOINT  SPACING 

ALSO  FIELDED  WITH  JOINTS  ORIENTED  AT 
0°  AND  90'  TO  W  P 


properties  in  the  continuum  approach  is  usually  hard  to  justify, 
and  this  approach  precludes  discontinuities  in  displacements 
across  joints,  which  may  be  a  critical  omission.  The  properties 
of  joints  in  the  explicit  approach  are  hard  to  measure,  and  the 
arrangement  of  joints  in  a  real  field  around  a  real  tunnel  can 
only  be  estimated. 

1.1  OBJECTIVE  AND  APPROACH 

Tin..  >b.uctive  of  this  study  is  to  examine  the  influence  of 
t he  joints  on  tunnel  closure  and  damage  using  numerical  techniques 
which  explicitly  treat  the  idealized  joints  in  the  SRI  structure. 

The  approach  involves  three  related  aspects: 

1 .  Joi n t  Ac t i vat ion . 

An  analytical  closed  form  solution  is  derived  for  the 
initiation  or  activation  of  relative  motion  across  a 
joint.  The  joint  activation  depends  on  (a)  the  rela¬ 
tive  orientation  of  the  joint  and  the  tunnel  surface 
intersected  by  the  joint,  { b )  the  magnitude  and  propa¬ 
gation  direction  of  the  stress  wave  encompassing  the 
tunnel.,  and  (c)  the  shear  stress  on  the  joint  which  must 
exceed  a  critical  value  which  depends  on  the  normal 
stress  across  the  joint  as  well  as  on  surface  roughness. 

tn  Section  J,  the  joint  activation  solution  is  derived 
and  applied  to  the  SRI  small-scale  models  to  estimate 
(a)  the  areas  on  the  tunnel  surface  where  joint 
activation  is  likely  to  occur,  (b)  the  stress  levels 
which  are  necessary  to  activate  joints  of  various 
or  1  cut a 1 1 ons ,  and  (c)  the  sensitivity  of  joint 
a  c  t..  i  v a  t  i o n  to  ma  t e  i  i  a  1  p  r  o }  >e  rti.es. 


2 .  Dynamic  Tunnel  Response  Cases . 


In  the  SRI  small-scale  models  with  tunnel  diameters  of 
3" ,  the  transit  time  across  this  characteristic  dimen¬ 
sion  is  only  22  usee  in  the  16A  rock  simulant  material. 
Tims,  since  the  rise  times  to  peak  stresses  of  0.8  kbar 
in  the  DIABLO  HAWK  event  are  A5-10  msec,  these  small- 
scale  structures  are  actually  loaded  (and  unloaded)  in 
a  quasi-static  fashion.  Full-scale  tunnels  may  be 
loaded  more  dynamically,  therefore  the  effects  of 
dynamic  loading  on  small-scale  tunnels  buried  in  jointed 
media  are  analyzed  using  six  plane  strain  computer  code 
simulations  with  dynamic  loading  times  of  100-400  s ec . 
These  cases  also  vary  the  assumed  joint  and  tensile 
strength  properties  of  the  rock  simulant.  The  WAVE-L 
Lag rang ian  code  is  used;  the  basic  formulation  is 
similar  to  the  HEMP  code  (Ref.  3).  The  numerical 
solutions  and  material  properties  are  described  in 
Section  3. 

3.  Ouasi-S tutic  Tunnel  Response  Prediction  For  a  Small-Scale 
Structure  Loaded  to  0.8  kbar  During  the  DIABLO  HAWK 
I  ',  ven  t . 

A  numerical  simulation  prediction  of  tunnel  closure  and 
damage  is  made  and  compared  with  a  fielded  small-scale 
API  r;oael  ( i  he  lias  a  3"  tunnel  diameter  with  q" 


An  quasi -st  m  is  solution  technique  and  the  caiculat  ronal- 
■  ■  r  l  mo  nt.  a  1  cops  at  isons  a  re  described  in  Section  4. 


ill  sons  are  described 


1.3  SUMMARY  AND  CONCLUSIONS 
1.3.1  Joint  Activation 


Although  the  response  of  a  tunnel  in  a  jointed  medium  must 

be  solved  numerically,  the  conditions  that  will  initiate  joint 

slipping  can  be  analyzed  in  closed  form  if  we  assume  an  elastic 

media  with  joints,  and  a  far  field  stress  state  of  uniaxial 

strain  with  a  stress  loading  of  o, 

load 

Figure  2  shows  the  loading  stress  (o,  J  required  to  activate 
a  joint  oriented  at  a  =  0,  4 45°,  and  90°  as  a  function  of  the 

angle  (3)  wnere  the  joint  intersects  the  tunnel  surface.  Thus 
a  joint  oriented  at  i  =  45°  will  require  a  loading  stress  of 
.06  kbar  (  900  psi)  to  become  activated  if  this  joint 
the  tunnel  surface  at  3  =  100°  (see  the  following  sketch  and 
Figure  2) . 


N 

N 

\ 


(Un  lined  tunnel  approximat ion  is  valid 
because  r\  ;n*.  :  involves  only 

small  strains  and  thus  the  backpacking 
in  the  tunnel  will  not  contribute  a 
significant  stress  to  the  tunnel  surface.) 


\  Symmetrical  Joint 

Behavior  for  3  +  180 

\ 
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Note  that  for  each  joint  orientation  (•*),  the  shaded  region 
on  Figure  2  shows  the  angles  (. )  on  the  tunnel  surface  whore  joint 
activation  is  possible  as  a  function  of  loading  stress  (  ^oac j)  • 

Note  that  symmetry  implies  equivalent  tunnel  behavior  at  angles 
of  .  +  180°.  Joints  will  not  slip  outside  of  these  shaded  regions. 

The  tunnel  angles,  .  ,  associated  with  absolute  minimums  in  load¬ 
ing  stress  ('  ,)  necessarv  for  joint  activation  are  indicated 

J  load 

in  Table  1.  However,  note  that  the  activation  stress  curves 
(Figure  2)  are  relatively  flat  bottomed;  and  therefore,  joint 
activation  can  occur  over  a  fairly  wide  range  in  tunnel  angles 
with  only  a  small  increase  in  loading. 

Table  1.  Tunnel  Angles  (r)  for  Minimum  Loading 

Stress  lor  Various  Joint  Orientations  (  t ) . 


Joint 

Orientation  I 

Tunnel 

ij 

Anglos  for  Minimum 
oading  Stress* 

i 

a eg  roe 

1  * 

l  degree 

kbar 

load 

psi 

0 

i 

6  7 ,  113 

.  068 

986 

4b  j 

103 

.059 

8  5b 

1 

41,  139 

.  1.62 

2350 

....  .1 

7b,  103 

.059 

8  5b 

*  Representative  constitutive  properties  are  used 
(Case  pl)  :  Joint  Cohesion  (io)=.04  kbar,  Friction 
Angie  (:)  31°,  Material  Poisson's  Ratio  (v)=.25. 

Figure  >  shows  the  areas  on  the  tunnel  surface  where  joint 
ict.  lvation  (tor  a  single  joint)  could  occur  for  a  .  -  .5  kbar. 
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All  the  joint  orientations  in  the  SRI  experiment  are  shown  for  a 
representative  set  of  joint  and  earth  media  properties.  (Varia¬ 
tions  in  the  joint  and  earth  media  constitutive  properties  are 
illustrated  in  Section  2.) 

1.3.2  Dynamic  Response  Cases 

Six  WAVE-L  numerical  simulations  (Table  3)  were  performed 
involving  3.00-400  nsec  stress  rise  times  to  a  peak  loading  stress 
of  0.S  kbar  on  the  SRI  3M  diameter  models.  The  wave  transit  time 
across  the  3"  tunnel  opening  is  22  |isec,  therefore,  the  dynamic 
tunnel  response  characteristics  could  be  investigated. 

Tile  tensile  stress  history  near  the  tunnel  crown  is  strongly 
dependent  on  i he  stress  wave  rise  time  for  times  corresponding 
to  a  t ew  transit  times  across  the  tunnel  diameter.  A  stress 
wave  rise  t  i  ;:k  oi  100  jisec  (-5  transit  times)  caused  a  peak  ten- 
si  U  .it  i  is'..;  almost  ini)  bars.  In  the  quasi-static  loading  case 

wiineat  i  ;it  l  Leinn,  no  tensile  stress  developed  near  the  crown. 

e.-u.  i  e  u  I'S.iti't  et  the  material  near  the  tunnel  surface  is 
uise  t<  i  t’i  ;  ?  ;u.  t  <  Lit  i  vo  displacements*  across  the  joints. 

Win  n  ..i  i;  ;  ,n  :  ■  *cai  n  across  joints,  the  plates  of  rock  separated 
;  ;  t  ..v  . :  ■  i : :  *  . :  can  act  as  unsupported  beams  due  to  the  presence 

•  u  t  ne  an  i  in^.a  tunnel  .  7a»  result  inn  beam  bonding  can  load  to 

*  i  I  ■  .  c  .  a  i  I  r.icture. 

1  .  .  i  n,ia:.n-St  atic  Tunnel  Response  Predict  ion 

Ine  ;ua.u -static  solution  technique  is  described  in  Section  4. 
bamorical  siinalat  ion  j  u  < \i  i  c  t  i  ons  using  this  technique  were  made 
t  n  a  1  i  o  l  at.  d  sinn  1  1  -  sea  1  e  (  \u  diameter  tunnel)  SRI  model.  The 

*Yii'.  ■  h’l.iliw  ■ Y  i  sq  >  1  ucemeii  t.  or  si  lppiiui  across  a  joint,  is  measured 
by  tun:  current  4i  stalin'  between  t.wo  points  which  were  initially 
touch  inn  each  other  on  opposite  laces  of  the  joint. 


calculation  predictions  were  made  prior  to  excavation  and  examin¬ 
ation  of  the  model. 

Figure  4  shows  a  comparison  of  the  predicted  tunnel  closure 
(as  a  %  of  the  diameter)  versus  tunnel  angle,  .  The  agreement 
is  reasonable  except  near  the  crown-invert  ( 3  =  0°  and  180°) . 

The  maximum  predicted  closure  is  1.2?.  at  o  =  30°.  The  experi¬ 
mental  data  show  a  maximum  1.6%  closure  at  the  crown-invert 
compared  to  a  0.5%  predicted  value  at  the  crown-invert.  This 
difference  may  be  due  to  the  moderate  sliding  on  the  joints  and 
associated  tensile  failure  and  relaxation  which  occurs  in  the 
calculation,  but  not  in  the  experiment. 

In  the  post-shot  experimental  model,  there  was  deformation 
of  a  thin  (.003  inch)  stainless  steel  tunnel  liner  which  surrounds 
the  foam  filler.  This  liner  may  have  inhibited  sliding  along  the 
joints.  The  observed  liner  deformation  is  associated  with  the 
joints  which  are  predicted  to  slide  the  most;  thus,  some  "block 
motion"  appears  to  have  occurred. 

It  is  concluded  that  tunnel  closure  in  a  well-charac terized 
jointed  media  could  be  predicted  using  existing  computer  code 
techniques  involving  explicit  descriptions  of  joint  locations, 
orientations,  and  friction  properties.  Also,  because  of  the 
important  role  joint  (and  fault)  motion  plays  in  determining 
damage  to  tunnel  liners  and  nearby  earth  material,  it  is  felt 
that  continuum  approaches  to  predicting  tunnel  closure  and  damage 
might  not  be  as  useful. 
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SECTION  2 


JOINT  ACTIVATION  ANALYSIS 


The  necessary 
ping  depend  on  the 


conditions  for  a  joint  or  block  to  begin  slip- 
relative  .on*,  k: ;Y*»:  of  the  joint,  the  tunnel 


surface,  and  the  stress  wave,  as  indicated  in  the  following 
sketch : 


t 

\ 

\ 


1  Oil  U 


s-  is  the  joint  orienta¬ 
tion  angle 

r  is  the  tunnel  angle 
(angle  where  joint 
in  terser ts  tunnel) 

.  is  the  angle  between 
t lie  joint  and  the 
normal  to  the  tunnel 
sur  1  are  ,  .  =  ;  .■ ■  -  >  \ 


Joint  slippage  will  occur  when  the  shear  stress  (i)  on  the  joint 
exceeds  a  critical  value,  i  .  This  value  defends  on  the  normal 
stress  (  ■  )  across  the  joint  as  well  as  on  the  surface  roughness. 
Prior  to  joint  activation,  the  earth  media  around  the  tunnel  is 
assumed  to  respond  elastically  and  statically  to  a  uniaxial  strain 
.load  from  a  (guasi-sta tic )  pLanar  wave. 


1  b 


i* 


The  key  characteristics  of  interest,  with  respect  to  Joint 
no  *  ;  :v;  :  :  on  are  : 

•  the  areas  on  the  tunnei  surface  where  joint 
activation  is  likely  to  occur, 

•  the  stress  levels  which  are  necessary  to 
activate  joints  of  various  orientations, 

•  the  sensitivity  of  joint  activation  to  material 
properties . 

Tne  determination  of  these  characteristics  involves  the 
foil owi ng  e  1  erne nts  : 

I.  The  a  :  a  :  .vo-o>i :  in  /  a  •  /*.  .*;•  field  on  the  unlined 
tunnel  surface  (this  stress  field  is  shown  in 
Figure  5  as  a  ratio  with  the  loading  stress, 

1*e"  0 /oload) 

2  .  The  • :  n  :  - .  n :  ■;  :  .<  I  i  .v  no.J  <  ( l  .  e  .  ,  we  assume 

i  =  i  +  (tan  1 ) a  where  i  and  o  are 
max  o  n  o  1 

experimentally  determined  parameters) . 

J.  The  ;7,  /  (nj  , 

:  '  -  ‘  as  a  function  of  joint  orientation 

(  0  ,  position  on  the  tunnel  (2) /  and  material 
properties . 

The  activation  stress  load  relationship  (derived  in 
Appendix  A)  is 


a 

load 


(2  i  ^cosy  )  /  (sm  (2:,./-  ;  )  -  sin  1) 


2  -  11  +  2 cos  ( 2^ )  ] 

[1~\) 


(1) 


where 


V  =  ja-^|  (see  sketen  on  previous  page) 

L  and  y  are  joint  constitutive  parameters 

is  the  Poisson's  Ratio  for  the  earth  material 
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Note  that  Equation  1  shows  that  t  f  and  v  are  the  con¬ 
stitutive  properties  of  interest  for  joint  activation.  Table  2 
lists  the  numerical  values  of  these  properties  which  were  used 
to  calculate  joint  activation  relationships. 

Table  2.  Joint  and  Earth  Media  Constitutive 
Properties 


[  .  i 

Joint 

i 

Earth  Media 

Case 

Cohesion 

*  T 

'  o 

Friction  Angle* 

Poisson  1 s  Ratio 

kbar 

psi 

4)  (degrees) 

V 

pi 

.04 

— 

580 

33 

.25 

P2 

.  08 

1160 

33 

.25 

P  3 

.02 

290 

33 

.25 

P4 

.04 

580 

0 

.25 

P5 

i 

.04 

1 _ 

580 

i 

33 

.40 

*Tne  joint  cohesion  (i  )  and  friction  angle  (!)  are  related  to 

the  unconfined  compressive  strength  (o  )  of  the  earth  media  by 

2  :  cos  I 
o 


Figure  2  and  Figures  o  to  9  show  the  variations  in  '  , 

3  3  load 

versus  tunnel  angle  ( ;* )  for  the  joint  and  earth  media  property 
variations  indicated  in  Table  2.  Figure  3  and  Figures  10  to  13 
show  the  corresponding  regions  on  the  tunnel  surface  where  joint 
activation  is  possible  for  o.  3 


.5  kbar. 
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r  i  on  tat  ions 
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Comparisons  of  the  results  from  the  constitutive  property 
cases  Pi  to  P5  leads  to  the  following  conclusions: 

1.  Increasing  (or  decreasing)  the  cohesion,  i  , 
increases  (or  decreases)  the  minimun  loading 
stress  for  joint  activation  while  leaving 
unchanged  che  tunnel  angles  (8)  where  joint 
activation  is  most  likely. 

2.  Reducing  the  fiction  angle  to  vf?  =  0,  thereby 
simulating  a  von  Mises  limiting  shear  strength 
of  t  ,  causes  an  expansion  of  the  possible  joint 
activation  region. 

3.  Varying  Poisson* s  Ratio  for  the  earth  media  from 
v  =  .25  to  v  =  .40  causes  an  expansion  of  the 
possible  joint  activation  region. 
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SECTION  3 


DYNAMIC  RESPONSE  CASES 

In  the  DIABLO  HAWK  experiment,  the  stress  wave  time  history 
involves  rise  times  to  a  peak  stress  of  0.8  kbars  in  b-10  msec. 
In  Section  4,  quasi-static  techniques  are  used  to  treat  this  rise 


Time  Alter  Wave  Arrival  at  O.N  kbar  Stress  Ran.ue 


time,  which  is  quite  loruj  with  respect  to  tunnel  response.  In 
Section  3.1,  material  properties  for  all  the  cases  are  described. 
In  Section  3.2,  six  numerical  calculations  involving  much  shorter 
rise  times  of  between  100  and  400  Msec  are  described  and  compared 


3.1  MATERIAL  PROPERTIES 


10 A  Rock  Simulant 


In  the  small-scale  experiments,  the  jointed  media  is  16A 
rock  simulant-  The  elas tic-plas tic  properties  are  modeled  with 
the  following  parameters: 


Bulk  Modulus  K 

Shear  Modulus  G 

Poisson’s  Ratio  v 

Density  a 

o 

Yield  Strength  Y 


- 1 3 1  Mba  r 
.086  Mbar 
.23 

2.0  gm/cm3 

.22  kbar  4*  1.2P,  where  P 
is  pressure.  A  Prandtl-Reuss 
non-associa ted  flow  rule  is 
assumed • 


Tensile  failure  is  modeled  in  most  of  the  solutions  using  a 
tensile  failure  model  which  permits  oriented  cracks  to  develop 
(Ref.  b ) .  A  300  psi  =  34.5  bar  tensile  strength  was  measured 
m  a  simple  tensile  test  (Ref.  7)  on  a  16A  rock  simulant  sample. 
Most  of  our  solutions  use  this  value;  in  one  of  the  cases,  the 
tensile  strength  is  reduced  to  6-9  bars  (100  psi)  as  a  sensitivity 
study . 

In  two  of  the  cases,  frictional  forces  along  the  sliding 
joints  are  included.  In  these  situations  the  joint  constitutive 
mode 1  i s  s pec i f i ed  by : 


'  max 


+  ; ; .  i 

o  n 


;  +  (tan; ) 

o 


w  he  r  e 


-  maximum  allowable  shear  stress  on  joint 

ma  x 
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a  =  normal  stress  on  joint 
n  J 

r  =  the  joint  cohesion 
o 

U  =  (tanj>)  is  the  coefficient  of  friction  and  $ 
is  the  angle  of  friction 

The  experimental  values  (Ref.  7)  are 

i  =10  psi 
o 

-  30° 

The  Drucker-Prager  plastic  yield  surface  assumed  for  the  16A 
rock  simulant  material  does  not  include  shear  strength  softening 
(i.e.,  reduction  of  the  yield  strength)  due  to  plastic  shear 
deformation  or  tensile  cracking.  This  may  be  significant  in  this 
application  because  the  uniaxial  strain  loading  path  for  this 
material  is  close  to  the  yield  surface.  therefore,  small  changes 
in  the  failure  surface  may  lead  to  large  changes  m  the  plastic 
response  oi  the  IGA  rock  simulant. 

Foam  Filler 

Trie  tunnel  is  filled  with  a  low  density  elastic-plastic  foam 
with  the  following  properties: 

Bulk  Modulus  K  =  1.75  kbar 


Shear  Modulus 
Density 

Von  Mises  Yield 
Strength 


G  =  2.5  kbar 
o  =  .05  gm/cirr 

Y  =  500  psi  =  0.34  kbar 
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3.2  CASE  COMPARISONS  FOR  THE  0.8  KBAR  DYNAMICALLY  LOADED  TUNNELS 


Table  3  lists  the  joint  geometry,  tensile  limit,  rise  time, 
and  material  model  for  the  six  dynamic  loading  cases.  In  all 
cases  the  tunnel  is  filled  with  low  density  foam  with  properties 
described  in  Section  3.1. 

The  loading  is  a  plane  wave  of  0.8  kbar  magnitude.  The 
stress  increases  linearly  from  0  to  0.8  kbar  over  a  characteris¬ 
tic  rise  time  t  ,  (100  psec  in  most  of  these  cases)  .  A  constant 

pressure  of  0.8  bar  is  maintained  after  t  .  The  far-field  con¬ 
dition  of  the  pulse  (i.e.,  prior  to  interaction  with  the  tunnel) 
is  uniaxial  strain. 

3.2.1  Effect  of  Loading  Rise  Time 

In  the  first  three  cases,  the  rock  media  is  modeled  as  an 
elas t lc-plas tic  media  with  no  joints  or  tensile  failure.  Cases 
1,  2,  and  3  differ  in  their  loading  rise  times  (t  =  100,  200, 
and  400  psec)  .  Cases  4,  5,  and  6  have  joints  which  are  perpe: 

dicular  to  the  wave  vector  as  indicated  in  Table  3;  in  these 
cases,  the  joint  properties  and  tensile  strength  of  the  media 
are  varied. 

Figure  14  shows  particle  velocity  (upper  part  of  figure) 
and  principal  stress  (lower  part)  fields  at  62  Msec  for  Case  1. 
Note  that  the  velocity  and  stress  fields  have  a  plane  of  symmetry 
which  separates  the  two  fields.  The  scale  bars  for  the  velocity 
vectors  and  principal  stresses  are  shown  on  this  figure.  The 
following  sketch  defines  the  plottinq  convention  for  the  principal 
s tress : 


Table  3.  The  Dynamic  Loading  Cases 


Cast: 

Joints 

Tensile  Limit 

Rise  Time  T 

i 

1 

None 

None 

100  Usee 

2 

None 

None 

200  usee 

3 

None 

None 

400  usee 

j 

1  ^ 

Frictionless 

Xto  Wave  Vector 

i 

100  psi  (6.9  bars) 

100  psec 

i 

Frictionless 

Xto  Wave  Vector 

300  psi  (34.5  bars) 

100  usee 

1 - - 

i  6 

With  Friction 
('max  =  10  psi  +  .6  on) 

X  to  Wave  Vector 

500  psi  (34.5  bars) 

100  usee 

St  res 
Wa  ve 


Cases  4  %  3  ,  f) 


critical  tensile 
in  t he  i nd i c a ted 
direction ,  '  >  . 


fa i 1 ure 
p  r i n c i pal  stress 


Plane  which  experiences  tensile 
f a i lure . 


Crack  opens  and  relieves  the 
tensile  stress  such  that  no  lorces 
or  t  rue t  i o n s  e x i s t  on  the  crack 
surfaces;  i.e.,  zero  normal  and 
shear  stress  on  1  he  crack  surface. 


In  the  directions  perpendicular 
to  the  c  r  ac  k  surf  a ce  ,  the  material 
can  maintain  tensile  stresses 
until  the  critical  tensile  failure? 
criterion  is  exceeded  and  a  new 
crack  is  formed. 


P 1  ot  t.  i  nq  Con  veil  t.  ion  : 


le  crack  in  x-y  plane, 
direction  of  the  crack 
ss  cams  i  no  the  crack)  . 

,  > !  -  pi  a  in:  crack  ;  i  .  o .  , 

■  c  -pi  am '  . i nd  x-y  piano 

x-y  plane  cracks;  only 
i  t  t  ed  . 


The  line  is  plotted  in 
(normal  to  the  principal 


caused  the  crack. 


ci  ack  . 


ou*-of-p]  me  tension  is 


1  ftely  !  i act  a i e » 1  -  no  tension  permitted. 


1*0 


lensilf  Con  ks  in  WAV1-1,  C.nlt'. 


I  to  ot  run  l 


In  Cases  4  and  5  the 


of  the  solution  to  varia¬ 


tions  in  the  tensile  strength  of  the  rock  medium  is  examined. 

In  Case  4,  the  IbA  rock  simulant  has  a  tensile  strength  of 
b.9  bars  (=  10  0  psi)  ;  in  Case  5,  a  strength  of  34.5  bars 

500  psi)  is  used.  The  effect  on  the  tensile  crack  formation 
is  shown  in  Figures  17  and  18  which  show  these  two  cases  at  two 
different  times. 

Case  0  implements  all  the  material  and  joint  properties  used 
in  thd  quusi-s tutic  DIABLO  HAWK  simulation  (Section  4),  including 
friction  along  the  joints.  The  presence  of  the  friction  is  the 
only  difference  between  Case  5  and  6.  The  final  crack  pattern  in 
Case  6  is  compared  to  Case  5  in  Figure  19.  i-'i*  '  '  ■>.  r  ;  •* 

Tensile  cracking  of  the  material,  near  the  tunnel  surface  is 
closely  related  to  the  relative  slide  displacements  across  the 
joints.  When  slipping  occurs  across  joints,  tile  pintes  o‘  rock 
separated  by  tile  joints  can  act  a s  unsupported  beams  due  to  the 
i^resence  of  uie  unlined  tunnel.  Tile  resulting  beam  bonding  leads 
to  tensile  stresses  ana  fracture. 


I'onsiU-  Slivn>\lh  -  t>.tJ  h.irs  -  100  j>s  i 


Tensile  Strength  =  6.9  bars  =  100  psi 


BiiwifllBHBSSSUaagBySgggllBligUilififiB 

^■■■■■■■■BS5S«ssi!SsyyssBB«8S!SSSB5SS5S 

S&SSSBSSSBifi5^2SS5s»Si2ss5&BiBa«BUBB 


hrHHHI 

igaag^isIliBgBBHBii 


DBQDI 


EiPiHli 


X  Cn 


Tensile  Strength  =  34.5  bars  =  500  psi 


-7  -0  -5  -4.  -'3  -2  -1  o  1  2  3  *i  5  6  7 

X  CH  I2e-5USEC 

iilM,r>.  iri.  tensile  Craes  Coni  igurat  ion  at  t  =  1-7  ,see  lor  lens  i  le  Strength 
,,|  too  psi  and  300  psi  (Cases  4  and  3;  O.S  khar  =  1  I • 6  ksi 
Sircss  Loading  in  100  i.soc)  . 


.  ’)  b.  U's 


'jOO  psi 


SECTION  4 


DIABLO  HAWK  SMALL-SCALE  STRUCTURE  PREDICTION  CASE 

This  numerical  simulation  was  intended  to  match  a  fielded 
small-scale  SRI  model.  The  calculation  results  were  obtained 
■  to  the  excavation  of  the  model.  The  particular  model  is 

designated  as  S-Dl-2-5  by  SRI. 

4.1  PROBLEM  CONFIGURATION 

The  tunnel  diameter  is  3"  .  Joint  spacing  is  nominally  '2" 

(0.47")  .  The  boundaries  of  the  numerical  grid  are  at  9.4  5"  from 

tile  tunnel  center  ( 0 . 3  times  the  tunnel  radius)  ,  which  is  suf¬ 
ficient  to  minimize  edge  el:  foots  in  the  region  of  interest.  The 

]  o  i  n  t  go  ome  t  r  y  is  s!  iown  w  i  t  ii  the  initial  comp  u  t  a  Lion  a  1  grid  in 
Figure  20. 

Material  properties  of  the  IGA  rock  simulant  and  foam  filler 
are  described  in  Section  3.1.  The  frictional  properties  ul  the 
joints  are  also  discussed  in  that  section. 

The  tunnel  and  jointed  media  are  subjected  to  a  planar  stress 
wave  of  0.8  kbur  peak  amplitude.  The  rise  Lime  is  5- 1 0  ns<-c 
and  decay  I'M.)  rnsi.c.  Tnose  times  are  very  lone  compared  to 
the  2  2  ..sec  sound  wave  transit  time  across  the  tunnel  diameter 
in  the  IGA  material.  The  actual  tunnel,  therefore,  is  in  a 
" guas l-s tat ic"  state. 

4.2  QUASI “STATIC  SOLUTION  TECHNIQUE 

The  juasi-static  solution  is  found  as  indicated  on  Figure  21: 

1.  Obtain  the  eLastic  solution  to  a  tunnel  loaded  with 
tile  0.8  kbar  (fnr-tield  uniaxial  strain)  stresses. 


4  2 


Figure  21.  Successive  Steps  in  Quasi-Static  Numerical  Simulation  Technique. 


2.  Using  tne  elastic  solution  as  initial  conditions, 
activate  the  elas tic-plas t ic  material  model  and 
calculate  the  dynamics  until  a  steady  state 
solution  is  obtained. 

3.  Using  this  elastic-plastic  solution  as  initial 
conditions,  a  calculation  is  performed  with  the 
slide  lines  introduced  to  model  the  joints  and 
permit  block  motion  and  deformation.  Also,  the 
material  may  fail  in  tension  if  the  tensile 
stresses  exceed  500  psi  in  any  principal 
direction.  The  dynamic  calculation  is  continued 
until  tensile  failure,  block  motion  and  tunnel 
deformation  have  stopped,  and  a  steady-state 
solution  is  thereby  obtained. 

4.3  JOINT  SLIDING 

Peak  slide  displacements  of  1Z  of  the  tunnel  diameter  are 
predicted  in  the  quasi-static  solution.  In  a  full  sized  tunnel, 
with  a  diameter  of  10  feet,  this  corresponds  to  1.2  inches  of 
sliding. 

Figure  22  shows  the  slide  displacement  along  the  joint  as  a 
function  of  distance  (scaled  by  tunnel  radius)  from  the  tunnel 
axis.  The  sliding  is  a  maximum  at  about  one  radius  along  the 
joints  nearest  the  tunnel  crown  and  invert  (Joints  3  and  9  in 
Figure  22) . 

4.4  TENS  I  LI-.;  CRACKING 

The  tensile  cracking  is  closely  related  to  the  relative  slide 
displacements  across  the  joints.  When  slipping  occurs  across  a 
joint,  the  plates  of  rock  separated  by  the  joints  act  as  beams. 
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The  presence  of  the  unlined  tunnel  removes  localized  support, 
and  beam  bending  leads  to  tensile  stresses  and  fracture. 


The  crack  pattern  obtained  in  the  numerical  solution  is 
shown  in  Figure  23.  The  "beams”  terminating  in  the  opening 
show  significant  cracking  of  the  unsupported  cantilevered  ends. 

4.5  TUNNEL  CLOSURE  AND  COMPARISON  WITH  EXPERIMENTAL  DATA 

The  tunnel  closure  predictions  are  shown  on  Figure  24 .  The 
closure  along  the  vertical  axis  (i.e.,  crown-invert)  is  approxi¬ 
mately  0.5%.  The  maximum  closure  of  1.2%  is  located  at  0  =  30°. 

Figure  24  shows  a  comparison  of  the  actual  and  predicted 
tunnel  closure  (as  a  %  of  the  diameter)  versus  tunnel  angle.  The 
agreement  is  quite  good  except  near  the  crown-invert  (0  =  0°  and 
180°).  The  experimental  data  show  a  1.6%  closure  at  the  crown- 
invert  compared  to  a  0.5%  predicted  value.  This  difference  is 
probably  due  to  the  moderate  sliding  on  the  joints  and  associated 
tensile  failure  and  relaxation  which  occurs  in  the  calculations, 
but  not  in  the  experiment. 

In  the  quasi-static  calculation,  tensile  cracking  occurs 
after  the  joints  are  ’‘unlocked"  and  joint  sliding  begins  (see 
Section  4.4  and  Figure  22).  Tensile  failure  is  not  evident  in 
the  post-shot  experimental  model  (Figure  25) .  If  the  metal  casing 
(which  is  not  modeled  in  the  calculation)  reduced  the  sliding  on 
the  joints,  or  if  the  actual  joints  are  stronger  (higher  cohesion 
and/or  higher  coefficient  of  friction)  than  used  in  the  calcula¬ 
tion,  then  there  would  be  loss  block  motion  in  the  experiment, 
and  therefore  there  would  also  be  less  tensile  failure. 

Figure  25  shows  the  post-shot  experimental  model.  Note  the 
deformation  of  the  thin  (.003  inch)  metal  casing  (or  liner) 
surround i no  the*  foam.  This  casing  may  have  inhibited  sliding 
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alone  tile  joints.  .  >  •*?-  i  ■  -  *  ,r.-  .*  ,**>'  ■  ,«  ’  /  ,  >. 

•  ,  »  -  .  Thus,  “block  motion"  has 

occurred  as  predicted.  A  quantitative  description  from  the 
small- scale  model  of  the  amount  of  relative  displacement  across 
the  joints  is  iiiqhly  desirable  for  comparisons  with  the  WAVl>L 
code  predictions  of  Fiqure  22. 
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APPENDIX  A 

JOINT  ACTIVATION  CONDITIONS  NEAR  A  CYLINDRICAL  TUNNEL 


The  relationships  between  the  principal  stress  components 
(oi,  o?)f  joint  and  tunnel  angles  (*,  15),  and  constitutive  acti¬ 
vation  properties  will  be  obtained  using  a  Mohr  diagram  (Ref.  8). 
The  joint  constitutive  model  is  specified  by 


or 


where 


max 


t  + 

o 


PO 


n 


max 


1  +  (tan  1)  e 

o  n 


(A-  1 ) 


=  the  maximum  allowable  shear  stress  on  the  joint 
max  J 

=  the  normal  stress  on  the  joint 
n  J 

I  =  the  joint  cohesion 
o 

m  =  tan  ;  is  the  coefficient  of  friction  and  ;  is 
the  angle  of  friction. 


Figure  A-l  (page  57)  is  a  Mohr  diagram  which  shows  (1)  the 

shear  stress  and  normal  stress  -  for  the  joint  constitutive 

n  J 

model  (Equation  A-l),  and  (2)  the  possible  shear  and  normal  stress 
states  corresponding  to  the  principal  stresses  1  and  0  4-  .  The 
following  sketch  (page  54)  indicates  the  geometrical  interpretation 
of  the  Mohr  diagram  in  terms  of  the  stresses  on  a  joint  near  a 
tunnel  surface: 
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Doing  the  algebra  and  solving  for  o? ,  we  find, 


>  t  =  ,  /sin(2y-;)  +  sin  A  ^  ^  [  o  / _ sin  9 _ \  (A- 3) 

1  l  sin  ( 2^- 1)  -  sin  j.y  u  Ism  (2^--;;.)  -  sin  T 

Equation  A-3  relates  the  two  principal  stresses  such  that 

an  oriented  joint  (with  angle  v  between  the  tunnel  normal  and 

the  joint)  will  lead  to  joint  activation.  A  more  physical  and 
useful  relationship  for  a  plane  stress  wave  engulfing  a  tunnel 
can  be  obtained  on  a  :  ’  an n - . 

On  an  •>;.  /. >•;- ■ ./  tunnel  surface,  one  principal  stress  component 
vanishes;  i.e.,  g  =0.  Thus,  the  tangential  stress,  a2  ,  on  the 
tunnel  surface  for  joint  activation  is 


s  m 


\ 


2t  cos 
o 


sin  ( 2  ,  - )  -  sin 


sin(2v~v)  -  sin 


(A-4) 


When  a  plane  wave  engulfs  a  tunnel,  however,  the  tangential  sc 
varies  with  the  angle  on  the  tunnel  surface  as  illustrated  in  the 
following  sketch: 


1  oa<t 


Kel.it  ionship  lift  won  I  aiif'on  L  i  a  1  Stross, 

m3  Loading  Stross,  ■*  ,,  on  Ttinnol  Surlaoo 

1  oad 


In  the  sketch,  A(3)  is  the  amplification  factor  of  the 
°load  wave*  A  (3)  depends  on  the  angle  3  on  the  tunnel  surface 
and  on  the  elastic  Poisson* s  Ratio  (v)  of  the  media.  The  deri¬ 
vation  for  A(3)  can  be  obtained  by  superposition  of  elastic 
solutions  derived  in  Reference  8.  Figure  5  (in  body  of  report) 
shows  A (3)  for  v  =  .25  and  .4. 

Since  o  i  =  0  and  a2  =  A(3),  °ioad'  EcIuatlon  shows  that 


load 


_  o2 
A(3) 


2t  cos  (1) 
o 


[sin(2^-4>)  -  sin  ]  A  ( 3 ) 


(A- 


where 


A(3)  =  2  -  [1  +  2cos  ( 23 )  ]  1-^ 


This  is  the  activation  stress  load  relationship  (Equation  1) 
used  in  the  report. 


(compress i vo ) 


REFERENCES 


1.  S.  E.  Blouin,  "Data  Report,  DATEX  II,"  AFWL-TR-69- 1 4 9  , 

January  1970. 

2.  M.  Rosenblatt  and  K.  N.  Kreyenhag en,  "Free  Field  Motions 
in  Jointed  Media,"  DNA  2769F,  November  1971. 

3.  M.  Rosenblatt,  G.  E .  Eggum,  Y.  M.  Ito  and  K.  N.  Kreyenhagen, 
"Interactions  of  Stress  Waves  with  Deeply-Buried  Tunnels  in 
Geologies  with  Cracks  and  Discontinuities,"  DNA  4398F, 

June  1975. 

4.  B.  S.  Holmes  and  H.  E.  Lindberg,  "Small-Model  Hardened 
Structures  Experiments  for  Diablo  Hawk,"  SR177-27, 

January  1977. 

5.  M.  L.  Wilkins,  "Calculation  of  Elastic-Plastic  Flow," 
ECRL-7322,  Revision  1,  January  1969. 

6.  M.  Rosenblatt,  G.  E.  Eggum,  L.  A.  Do Angelo  and  K.  N. 
Kreyenhagen,  "Numerical  Investigation  of  Water  Drop  Erosion 
Mechanisms  in  Infrared-Transparent  Materials,"  AFML-TR- 76- 19 3 , 
December  1976. 

7.  D.  Stowe,  Personal  Communication,  U.S.  Army  Engineer  Water¬ 
ways  Experiment  Station,  Vicksburg,  MS. 

8.  S.  Timoshenko  and  J.  Goddier,  Theory  of  Elasticity,  Third 
Edition,  McGraw  Hill,  1970. 


58 


DEPARTMENT  OF  THE  AIR  FORCE  (Continued) 


DEPARTMENT  OF  ENERGY  CONTRACTORS  {Con tin uedj 


Air  Force  Weapons  Laboratory 

Los  Alamos 

National  Scientific  Lab 

Air  Force  Systems  Command 

ATTN 

B.  Killian 

ATTN:  NTE,  M.  Plainondon 

ATTN 

J.  Johnson 

ATTN:  NT*  l).  Payton 

ATTN 

MS  364 

ATTN:  SUL 

ATTN 

l. .  Germaine 

Air  University  Library 

Sandia  Laboratories 

Department  of  the  Air  Force 

L  i  v  e  rm  o  r  e  L  a  b  o  ra  t  o  r  y 

ATTN:  AUL/LSE 

ATTN 

Library  &  Security  Classificat 

As s i s ta n t  Chi e f  of  S ta f f 

Sandia  National  Lab 

Studies  A  Analyses 

ATTN 

L.  Hill 

De pa i' tme n t  of  t he  Air  f  o rc e 

ATI  N 

3141 

ATTN:  Af /SASM 

ATTN:  Af / SASM*  W.  Adams 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS 

Ballistic  Missile  Office 

Aerospace  Corp. 

Air  force  Systems  Command 

ATTN 

P.  Mathur 

ATTN:’  MNN 

ATTN 

Technical  Information  Services 

ATTN:  MNNH 

Aqbabian  Associates 

De p u ty  Chief  of  S ta  f f 

ATTN 

M.  Balachanda 

Operations  Plans  and  Readiness 

ATTN 

C.  Baqqe 

Department  of  the  Air  Force 

2  cy  ATTN 

M.  Aqbabian 

ATTN:  AfXODC 

Applied  Theory,  Inc. 

Deputy  Chief  of  Staff 

2  cy  ATTN 

J.  Trulio 

Research,  Development.  X  Acq 

Department  of  the  Air  force 

AVCO  Research  &  Systems  Group 

ATTN:  AFRDQ! 

ATTN 

Library  A830 

Foreign  Techno lon.v  Division 

BDM  Corp. 

Air  force  Systems  Command 

ATTN 

Corporate  Library 

ATTN:  NIIS  Library 

ATTN 

T.  Neighbors 

C  t  r  a  t.  e  <  |  i  l  Air  i  onm  a  n  d 

Booinq  Co. 

Department  of  the  Air  f  ori.e 

ATTN 

T.  Berg 

ATTN:  XPFS 

ATTN 

R.  Pyrdahl 

ATTN:  NRI-STINfO  Library 

ATTN 

M/S  4?/ 37 ,  K.  F riddel  1 

ATTN 

H.  Leistner 

OTHER  GOVERNMENT  AGiNCifS 

AT1N 

Aerospace  Library 

ATTN 

J.  Wooster 

Centra  1  Intelligence  Agency 

ATTN:  OSWR/NfD 

Cal ifornia 

Institute  of  Technology 

ATTN 

D.  Anderson 

Department  of  the  Interior 

Bureau  ot  Mines 

Cal ifornia 

Research  A  Technology,  Inc. 

ATTN:  Tech  l  ib 

ATTN 

1 ibrary 

ATTN 

S.  Schuster 

Department  fit  the  Interior 

ATTN 

K.  K reyen ha gen 

1 1 . .  Geol  Of)  i  ca  1  Survey 

4  cy  ATTN 

M.  Rosenblatt 

ATTN:  R.  Carroll 

4  cy  ATTN 

L.  De Angelo 

ATTN:  W.  Twenhofel 

Cal ifornia 

Research  A  Technology,  Inc. 

I )e [ -a r tme n t  fi f  the  In t e r  i o r 

ATTN 

D.  Orphal 

i.t.s.  Geological  Survey 

ATTN:  [).  Roddy 

University 

of  California 

ATTN 

N.  Cook 

Df  PAR IMF  NT  OF  ENERGY  CONTRACTORS 

ATTN 

R .  Goodman 

1  awrenre  L  i  vermore  Na  t.  i ona  1  l  ab 

Cal  span  Corp 

ATT  i :  'Vchniral  [nfo  Depr.  Library 

ATTN 

Library 

ATTN:  !. -PI,  0.  Oakley 

ATTN;  [,.  Smith 

University  of  Denver 

ATTN:  H .  Heard 

ATTN 

J.  Wisotski 

Oa  F  R  i  dqe  Na  t i ona 1  la  bo ra  t o  r  y 

EG&G  Wash  Analytical  Svcs  Ctr,  Inc. 

ATTN:  Central  Research  Library 

ATTN 

L ibrary 

60 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

DEPARTMENT 

OF  DEFENSE  CONTRACTORS  (Continued) 

Electromech  Sys  of  New  Mexico,  Inc 

R  &  D  Associates 

ATTN:  R.  Shunk 

ATTN 

D.  Rawson 

ATTN 

D.  Shrinivasa 

University  of  New  Mexico 

ATTN 

R.  Port 

Civil  Enqineerinq  Rsch  Fac 

ATTN 

J.  Lewis 

ATTN:  N.  Baum 

ATTN 

Technical  Information  Center 

ATTN 

P.  Haas 

Foster-Miller  Associates,  Inc 

ATTN:  J.  H amp son  for  E.  Foster 

Rand  Corp 

ATTN 

A.  laupa 

Franklin  Institute 

ATTN:  Z.  Zudans 

Sc ience  Appl ications,  Inc 

ATTN 

Technical  Library 

I l T  Research  Institute 

ATTN:  M.  Johnson 

Science  Applications,  Inc 

ATTN:  R.  Welch 

A T  T  N :  Doc  umen t  s  Lib ra r y 

ATTN 

Technical  Library 

Science  Applications,  Inc 

Institute  for  Defense  Analyses 

A 1  TN 

W.  Lay son 

ATTN:  Classified  Library 

Southwest  Research  Institute 

.  •! .  Wiqqins  iO,  Inc 

ATTN 

A.  Wenzel 

ATTN:  J.  Collins 

ATTN 

W.  Baker 

Kama n  AviDyno 

SRI  Internat ional 

ATtm;  Library 

ATTN 

G.  Abraham son 

ATTN 

H.  Lindberq 

Kaman  Sciences  .orp 

ATTN 

B.  Holmes 

A'*’  *  N  ;  l  i  I? rat* 

Systems,  Science  ft  Software,  Inc 

Kaman  TJMPC 

ATTN 

R.  Duff 

ATTN:  DA1*  I  At. 

ATTN 

C.  Arc hornbeam 

ATTN 

1.  ibrary 

!.  D.  -M  1  t  iwanqer  Cons  *  1  t  . 

ATTN 

D.  G  rifie 

A  T  ’N  :  ;  .  **a  1  t  iwanqer 

Terra  ’’’ok. 

Inc 

lOcFhe^d  Miss  )!(>•.  ft  '  ;*a<  e  Co,  Inc 

ATTN 

H.  Pratt 

ATrN:  Tec  nn  i r  a  1  I  n*  pmj  t  :  ;*n  l  enter 

ATTN 

L ibrary 

' 

Texas  A  ft  M  University  System 

Ma  s  sa  c  n  •  j  setts  Ins*  *  *  »'<  n  •  >  1  on  v 

ATTN 

A.  Rychlik 

ATTN:  e 

ATTN 

J.  Hand  in 

ve  r  ►'  1 1  f  ■  a  ->r  - ,  I  n*. 

TRW  Defense  ft  Space  Sys  Group 

A'TN:  ’ .  V»  • » •  *  i  *  f 

ATTN 

P.  Huff 

ATTN 

Technical  Information  Center 

1  :  /  )i  >;e  o *  V>w 

ATTN 

N.  Lipner 

TRW  Defense  ft  Space  Sys  Group 

N-.  hwe  s  torn  in  *  .v*''  i  *  * 

ATTN 

F .  Wonq 

ATTN 

P.  Dai 

Universal 

tnalvtics.  Inc 

V*,:  a.  i.;r;vP‘ 

ATTN 

['.  Field 

:  a<  i f  ;  .  - .  ierr  i  -esea c'  ►>  • « »* : ■ 

We idl inqer 

Associates,  Consulting  Engineers 

A '  "N  ;  ' ’ .  lormlov 

ATTN 

M.  Baron 

ArTN 

I .  Sandler 

'  \  f  i  (.a  Tec  nn >  ■ 1 .  q  v 

AT"N; 

We idl inqer 

Associates,  Consulting  Engineers 

ATTN 

.1 .  Isenberq 

;,atol  f  ner;ir  :  ses  ,  in- 

A'TN:  v. 

William  Ferret 

ATTN 

W .  Perret 

t  0 rri<i  *  i .  -fi,i  I  (: 

■  . 

*  i on  ;  ,» 1  ,i  hrary 

;  .  ,.1'jf  r 


61 


